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10  long-standing  problems  in  the  atmospheric  sciences  have  been  the  correct 
modeling  of  the  ion  chemistry  in  the  earth's  atmosphere  and  the  proper  deter¬ 
mination  of  the  ion  species  and  densities  through  in  situ  measurements.  Com¬ 
parison  between  experimental  data  and  simultaions  of  those  data  by  computer 
modeling  of  atmospheric  chemistry  is  a  means  of  validating  the  computer  code 
as  well  as  indicating  which  processes  are  in  need  of  further  study.  The 
DAIRCHEM  computer  code  is  used  here  to  simulate  data  taken  in  the  midlatitude  D 
region  during  quiet  conditions.  Comparison  between  the  total  positive  ion  - 
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ABSTRACT  (cont) 


ensity  profile  derived  from  rocket  measurements  and  the  one  computed  by  the 
code  shows  very  good  agreement  in  the  30  to  90  km  range,  with  the  exception 
that  the  simulated  ion  density  profile  is  of  somewhat  smaller  magnitude  than 
the  experimental  one  in  the  60  and  75  km  region.  Such  a  discrepancy  is  only 
partially  explained  by  the  uncertainties  in  the  ionization  rate  of  NO  due  to 
Lyman  alpha  radiation.  Comparison  between  the  measured  and  the  computed 
electron  density  profiles  shows  that  the  measured  profile  is  consistently  of 
a  smaller  magnitude  than  the  computer  profile  in  the  65  and  85  km  range. kWe 
interpret  this  discrepancy  as  a  deficiency  in  the  modeling  of  the  negatiPe"  ion 
chemistry.  Also,  this  deficiency  is  probably  the  main  cause  of  the  disparity 
between  the  total  positive  ion  density  profiles  in  the  corresponding  altitude 
range.  The  authors  feel  that  the  positive  ion  chemistry  is  reasonably  well 
understood  and  that  this  section  of  the  DAIRCHEM  code  may  be  considered 
validated.  The  negative  ion  chemistry  needs  to  be  studied  further.  Specifi 
cally,  alternate  electron  attachment/detachment  processes  should  be  considered, 
as  well  as  an  as-yet-undetermined ,  possibly  very  massive,  negative  species 
which  may  affect  the  ion  recombination  rates.  Successful  modeling  of  the  D 
region  under  naturally  occurring  conditions  is  one  means  of  validating  the 
computer  codes  which  are  used  as  input  for  Army  communications  systems  codes 
and  nuclear  weapons  effects  codes. 
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SUMMARY 


Experimental  data  from  rocket  measurements  taken  in  the  quiet,  daytime 
D  region  and  simulations  of  these  data  by  the  DAIRCHEM  computer  code  are 
compared.  The  purpose  Is  to  validate  the  computer  code  and  to  identify 
those  physical  processes  that  need  further  study. 

The  direct  comparison  of  simulated  and  experimental  data  is  shown  graph¬ 
ically.  It  Is  Important  to  note  that  no  data  have  been  normalized.  The 
total  positive  ion  density  profiles  show  good  agreement  over  the  30  to 
90  km  region  with  the  exception  of  the  60  to  75  km  span.  Possible  vari¬ 
ations  in  the  absorption  of  Lyman  alpha  radiation  could  account  for  part 
of  the  discrepancy,  although  it  is  felt  that  uncertainties  in  the  for¬ 
mation  of  negative  ions  and  the  total  negative  ion  and  electron  densities 
are  the  cause.  Qualitative  agreement  between  predictions  of  the  detailed 
positive  ion  composition  and  in  situ  mass  spectrometric  measurements  indi¬ 
cate  that  the  positive  Ion  chemistry  is  being  correctly  modeled.  The 
computed  dominant  positive  and  negative  Ions  are  listed  as  a  function  of 
altitude  In  a  table  which  serves  as  a  guide  when  comparing  these  results 
with  other  measurements  or  theoretical  calculations. 

Comparison  of  the  simulated  and  experimental  electron  density  profiles 
shows  that  the  experimentally  derived  values  consistently  lie  below  the 
predicted  ones  in  the  60  to  85  km  region.  Mass  spectrometer  measurements 
of  negative  ions  (currently  available  only  above  70  km)  show  little  agree¬ 
ment  with  current  predictions  In  the  70  to  90  km  range.  These  discrepan¬ 
cies  indicate  that  there  are  electron  attachment  mechanisms  and  one  or 
more  negative  species,  possibly  small  particulates,  which  are  not  being 
correctly  modeled. 
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INTRODUCTION 


The  Air  Chemistry  (AIRCHEM)  computer  code1  has  been  developed  to  simul¬ 
taneously  solve  a  large  number  of  stiff,  time-dependent,  differential 
equations2  and  thus  correctly  describe  the  ion  chemistry  associated  with 
the  nuclearly  disturbed  atmosphere.3  A  variation  of  this  code,  D  Region 
Air  Chemistry  (DAIRCHEM),  has  also  been  developed  to  model  the  chemistry 
of  the  ionospheric  region  under  both  quiet  and  naturally  disturbed  condi¬ 
tions.4  The  plan  is  to  compare  the  predictions  of  the  DAIRCHEM  code 
against  experimental  data.  Reasonable  agreement  between  the  measured 
data  and  simulations  of  that  data  by  the  DAIRCHEM  code  serves  to  vali¬ 
date  the  code,  whereas  disagreement  indicates  which  atmospheric  processes 
need  further  study.  Validation  of  the  DAIRCHEM  code  increases  the  level 
of  confidence  in  the  AIRCHEM  code's  predictions  for  the  nuclearly  dis¬ 
turbed  atmosphere  and  its  potential  effects  on  Army  systems. 

The  case  studied  here  is  the  daytime,  quiet,  midlatitude  D  region.  A 
series  of  coordinated  measurements  was  made  from  Wallops  Island  (37.8°  N, 
75.5°  W)  for  an  extended  period  of  time  which  covered  both  quiet  and  dis¬ 
turbed  days.  The  data  selected  for  comparison  were  taken  on  31  January 
1972,  a  quiet  day  with  no  anomalous  D  region  radio  absorption.  Two  rocket 
flights  at  12:30  and  13:05  local  time  (solar  zenith  angles  of  55°  and  58°, 


lE.  L.  Lortie,  M.  D.  Kregel ,  and  F.  E.  Niles,  1976,  "AIRCHEM:  A  Com¬ 
putational  Technique  for  Modeling  the  Chemistry  of  the  Atmosphere,"  BRL 
Report  1913  (AD  A030157) 

2T.  P.  Coffee,  J.  M.  Heimerl ,  and  M.  D.  Kregel,  "A  Numerical  Method  to 
Integrate  Large  Stiff  Systems  of  Ordinary  Differential  Equations,"  sub¬ 
mitted  to  BRL  for  publication 

3 (see  for  example)  J.  M.  Heimerl  and  F.  E.  Niles,  1978,  "BENCHMARK-76: 

Model  Computations  for  Disturbed  Atmospheric  Conditions  III.  Results 
for  Selected  Excitation  Parameters  at  60  km,"  BRL  Technical  Report  ARBRL- 
TR-02051  (and  references  therein). 

4D.  W.  Hoock  and  M.  G.  Heaps,  1978,  "DAIRCHEM:  A  Computer  Code  to  Model 
Ionization/Deionization  Processes  and  Chemistry  in  the  Middle  Atmo¬ 
sphere,"  Atmospheric  Sciences  Laboratory  Internal  Report,  White  Sands  Mlssil 
Range,  NM 


respectively)  provided  independent  experimental  measurements  of  the 
electron  and  total  positive  ion  densities. 5»6*7 


COMPARISONS  OF  EXPERIMENTAL  DATA  WITH  THE  DAIRCHEM  SIMULATIONS 

The  important  sources  of  ionization  for  the  D  region  are  X-rays,  Lyman 
alpha  radiation  (Lya)  and  galactic  cosmic  rays  (GCR).  The  ion-pair 
production  rates  (q)  used  in  the  computations8  are  shown  in  figure  1; 
the  sources  were  not  directly  monitored  during  the  rocket  measurements. 
Because  31  January  was  a  quiet  day,  the  X-ray  flux  was  not  large,  and 
the  direct  ionization  of  NO  by  Lya  was  the  largest  source  of  electrons 
throughout  the  D  region.  Due  to  the  rapid  attenuation  of  Lya  by  02, 

the  dominant  source  of  ionization  near  65  km  and  lower  becomes  the  GCR 
background.  The  source  of  ionization  does  affect  the  type  of  positive 

ion  which  is  produced.  During  quiet  periods,  N0+  is  the  most  abundant 
ion  initially  produced  throughout  the  D  region.  At  65  km  and  below,  the 

initially  produced  Ions  are  N^,  02,  N+,  and  0+,  in  that  order;9  but  due 

to  rapid  charge  exchange  and  ion-atom  interchange  processes,  the  effective 

ion  production  rate  Is  approximately  0.8  02  and  0.2  N0+  ions  per  ion  pair. 

The  effects  on  the  positive  ion  chemistry  will  be  discussed  in  the  next 
section. 


5L.  C.  Hale,  1974,  "Positive  Ions  in  the  Mesosphere,"  COSPAR  Methods  of 
Measurements  and  Results  of  Lower  Ionospheric  Structure,  K.  Rower,  ed. , 
219-235,  Akademle-Verlag,  Berlin.  See  also  J.  D.  Mitchell,  "An  Experi¬ 
mental  Investigation  of  Mesospheric  Ionization,"  Ionospheric  Research 
Scientific  Report  416,  Pennsylvania  State  University,  27  June  1973;  and 
T.  W.  Lai,  "Electron  Collection  Theory  for  a  D-Region  Subsonic  Blunt 
Electrostatic  Probe,"  Ionospheric  Research  Scientific  Report  424,  Penn¬ 
sylvania  State  University,  20  May  1974 

6L.  G.  Smith  and  K.  L.  Miller,  1978,  "The  Measurement  of  02  Number 
Density  by  Absorption  of  Lyman  Alpha,"  J  Geophys  Res,  79:1965-1968 

7E.  A.  Mechtly,  1974,  "Accuracy  of  Rocket  Measurements  of  Lower  Iono¬ 
spheric  Electron  Concentrations,"  Radio  Sci ,  9:373-378.  See  also  Progress 
Report  73-1,  p.  62ff,  Research  in  Aeronomy  October  1,  1972  -  March  31, 

1973,  edited  by  Belva  Edwards,  University  of  Illinois,  Urbana,  Illinois, 
1973 

8I.  Chldsey,  1978,  "SOURCE:  A  FORTRAN- IV  Subroutine  Package  that  Models 
Several  Naturally  Occurring  Energy  Inputs  to  the  Mesosphere,"  BRL  Technical 
Report,  ARBRL-TR-02093 

9F.  Gilmore,  as  quoted  in  table  6  in  B.  F.  Myers  and  M.  R.  Schoonover, 
"Electron  Energy  Degradation  In  the  Atmosphere:  Consequent  Species  and 
Energy  Densities,  Electron  Flux,  and  Radiation  Spectra,"  DNA  3513T, 

3  January  1975 
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ALTITUDE  (km) 


The  electron  and  total  positive  ion  density  profiles  computed  by  using 
ion  pair  production  rates  from  figure  1  are  shown  as  circles  in  figure 
2.  (Note  that  the  upper  scale  for  the  positive  ions  has  been  adjusted 
to  the  right.)  The  computed  electron  densities  and  total  positive  ion 
densities  are  essentially  equal  to  one  another  above  65  km.  Also  shown 
in  figure  2  for  comparison  are  the  total  positive  ion  density  profile 
(derived  from  blunt  probe  conductivity  data),  the  inferred  electron  den¬ 
sities  from  the  same  data,  and  the  measured  electron  density  profile  (com¬ 
bined  probe  and  propagation  experiment).  We  emphasize  that  the  profiles 
have  been  simply  overlaid;  no  data  have  been  normalized. 

Two  general  points  should  be  noted  when  the  profiles  are  compared.  First, 
the  total  positive  ion  density  profiles  agree  very  well  over  a  large  alti¬ 
tude  range,  with  the  exception  of  the  60  to  75  km  span  where  the  simulated 
total  positive  ion  density  profile  is  less  than  that  determined  from  the 
conductivity  measurements.  Second,  the  simulated  electron  density  profile 
is  consistently  greater  than  the  measured  and  inferred  profiles  in  the 
60  to  85  km  region.  (Measured  electron  densities  were  not  obtained  below 
60  km. ) 

Several  suggestions  may  be  initially  tendered  as  to  why  the  differences 
between  the  computed  and  experimental  profiles  exist.  The  smaller  simu¬ 
lated  total  positive  ion  densities  in  the  60  to  75  km  region  may  be 
partially  due  to  an  incorrect  determination  of  the  absorption  of  Lyman 
alpha  and  subsequent  ionization  of  NO.  That  is,  the  Lyman  alpha  may  have 
actually  penetrated  deeper  than  was  assumed  in  the  model  computation. 

This  could  happen  if  either  the  0 ^  column  density  were  less6  or  the  O2 

absorption  cross  section  were  smaller10  than  the  values  actually  used, 

-20  2 

i.e.,  1  x  10  cm.  Our  calculations  show  that  either  occurrence  improves 
the  agreement  between  the  total  positive  ion  density  profiles  at  70  and 
75  km.  (Unfortunately,  each  occurrence  also  increases  the  differences  in 
the  electron  density  profiles  in  the  same  range.)  However,  differences 
at  60  and  65  km  still  remain.  Another  possibility  that  may  explain  the 
differences  in  the  positive  ion  profiles  is  that  the  effective  rate  of 
recombination  may  be  too  large,  thus  reducing  the  computed  ion  densities. 
This  Is  related  to  the  fact  that  the  computer  code  predicts  more  elec¬ 
trons  than  the  measurements  indicate.  Finally,  this  tendency  toward 
larger  simulated  electron  densities  implies  that  there  may  be  one  or 
more  negative  ion  species  in  the  60  to  85  km  region  which  are  not  being 
properly  considered. 


6L.  G.  Smith  and  K.  L.  Miller,  1974,  "The  Measurement  of  O2  Number  Density 
by  Absorption  of  Lyman  Alpha,"  J  Geophys  Res,  79:1965-1968 

10J.  H.  Carver,  H.  P.  Gles,  T.  I.  Hobbs,  B.  R.  Lewis,  and  D.  G.  McCoy, 
1977,  "Temperature  Dependence  of  the  Molecular  Oxygen  Photoabsorption 
Cross  Section  Near  the  H  Lyman  Alpha  Line,"  J  Geophys  Res.  82:1955-1960 
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Figure  2.  Open  symbols  are  the  electron  density  profiles  (bottom  scale). 

Measurements  were  made  by  the  Pennsylvania  State  University  Group 
(ref  5)  and  the  University  of  Illinois  Group  (ref  7).  Solid 
symbols  are  the  total  positive  ion  profiles  (top  scale). 
Measurements  were  made  t.  the  Pennsylvania  State  University 
Group  (ref  5). 


11 


POSITIVE  ION  CHEMISTRY 


The  measurement11  more  than  a  decade  ago  of  water  cluster  ions  in  the  D 
region  substantially  revised  our  understanding  of  that  portion  of  the 
ionosphere.  Subsequent  laboratory  studies12’1"  showed  that  pathways  do 

exist  for  converting  ol  and  N0+  ions  to  water  cluster  ions.  The  essential 

+  ^  + 
points  of  the  02  reaction  path  are  the  three-body  clustering  of  02  with  02 

0*  +  02  +  M  -v  oj  +  M, 

a  rapid  switching  reaction  with  water 

o*  +  h2o  -  o2  (h2o)  +  o2, 

and  the  rapid  two-body  reactions 

02  (H20)  +  H20  -  H30+  (OH), 

H30+  (OH)  +  H20  H30+  (H20)  +  OH. 

Subsequent  hydration  to  form  H30+  (H20)n  may  also  occur.  While  this 
reaction  scheme  is  rapid  enough  to  explain  the  water  cluster  ion  concen¬ 
trations  where  02  is  the  major  precursor  ion,  such  as  below  65  km  where 

cosmic  rays  dominate  or  in  the  disturbed  D  region  when  ionization  by 
precipitating  particles  is  large,14  it  does  not  suffice  for  the  quiet 

time  D  region  where  N0+  is  the  major  precursor  ion.  Direct  sequential 

clustering  of  water  to  N0+  was  too  slow  to  produce  the  observed  water 


nR.  S.  Narcisi  and  A.  D.  Bailey,  1965,  "Mass  Spectrometric  Measurement 
of  Positive  Ions  at  Altitudes  from  64  to  112  Kilometers,"  J  Geophys  Res, 
70:3687-3700 

12F.  C.  Fehsenfeld  and  E.  E.  Ferguson,  1969,  "Origin  of  Water  Cluster  Ions 
in  the  D  Region,"  J  Geophys  Res,  74:2217-2222 

13F.  C.  Fehsenfeld,  M.  Mosesman,  and  E.  E.  Ferguson,  1971,  "Ion  Molecule 
Reactions  in  an  02  -  H20  System,"  J  Chem  Phys,  55:2115-2119 

14G.  C.  Reid,  1972,  "The  D  Region  During  PCA  Conditions,"  Magnetosphere 
Ionospheric  Interactions,  K.  Folkstad,  ed.,  39-46,  Universitets- 
forlaget,  Oslo 


cluster  densities.  More  rapid  switching  reactions  have  been  proposed15*16 
of  the  type 


N0+  +  N2  +  M  N0+  (N2)  +  M, 

N0+  (N2)  +  C02  -  N0+  (C02)  +  N2, 

N0+  (C02)  +  H20  -  N0+  (H20)  +  C02, 

with  a  similar  sequence  until  the  third  hydrate  is  produced.  The  third 
hydrate  then  undergoes  a  rapid  two-body  reaction  with  water 

no+  (h2o)3  +  h2o  -v  h3o+  (h2o)2  +  hno2, 

and  the  hydronium  ion  sequence  is  entered.17*18 

Repeated  rocket  measurements  over  a  wide  variety  of  latitudes  and  seasons 
have  shown  that  the  hydronium  ion  clusters  H30+  (H20)n  begin  to  rapidly 


15F.  E.  Niles  and  J.  M.  Heimerl,  1972,  "Association,  Switching  and  Re¬ 
arrangement  for  Positively  Charged  Cluster  Ions  in  the  Upper  Atmosphere, 
I:  Qualitative  Description,"  BRL  Report  1595 

160.  M.  Heimerl,  0.  A.  Vanderhoff,  L.  J.  Pucket,  G.  E.  Keller,  and  F.  E. 
Niles,  1972,  "Association,  Switching  and  Rearrangement  for  Positively 
Charged  Cluster  Ions,  II:  Applications  at  80  km,"  BRL  Report  1605 

17G.  C.  Reid,  1976,  "Ion  Chemistry  in  the  D  Region,"  Advances  in  Atomic 
and  Molecular  Physics,  12:375-414,  Academic  Press,  New  York 

10G.  C.  Reid,  1977,  "The  Production  of  Water-Cluster  Positive  Ions  in 
the  Quiet  Daytime  D  Region,"  Planet  Space  Sci,  25:275-290 


decrease  in  density  above  about  75  to  85  km.11*19"21  The  main  reason  is 
that  the  clustering  sequence  involves  several  three-body  reactions  with 
major  constituents  and  therefore  rapidly  decreases  in  efficiency  with 

increasing  altitude.  The  exact  shape  of  the  H.jO+  water  cluster  ion  pro¬ 
file  is  also  somewhat  a  function  of  the  temperature  and  water  vapor  pro¬ 
files.  Above  75  to  85  km,  one  might  expect  to  see  a  few  of  the  "early" 

clusters  in  the  N0+  and  0^  ion  hydration  chains. 

The  simulated  positive  ion  density  profiles  are  shown  in  figures  3  and  4 
for  ions  having,  respectively,  N0+  and  H30+  as  the  core  ions.  Also  shown 
is  the  predicted  total  positive  ion  density  profile.  As  expected,  the 
computer  code  predicts  that  N0+  and  its  first  two  water  clusters,  along 
with  the  intermediate  cluster  N0+  (C02),  dominate  the  upper  D  Region. 

There  is  good  qualitative  agreement  between  these-  predictions  and  several 
mass  spectrometric  measurements;19-22  nevertheless,  these  clusters  are 
weakly  bound  and  therefore  subject  to  collisional  breakup  from  the  in  situ 
rocket  measurement  techniques.  Figure  4  shows  the  computed  density  pro¬ 
files  for  the  H^0+  hydrates.  Above  65  km,  entry  to  the  hydronium  ion 
chain  is  through  NO  producing  H^O  (H20)2>  while  below  65  km,  entry  is 
through  0^  producing  H30+  (H20).  The  code  predicts  that  clustering  to 

heavier  hydrates  will  be  important  below  75  km.  Again,  the  comparison 
with  mass  spectrometric  measurements  gives  good  qualitative  agreement. 


nR.  S.  Narcisi  and  A.  D.  Bailey,  1965,  "Mass  Spectrometric  Measurements 
of  Positive  Ions  at  Altitudes  from  64  to  112  Kilometers,"  J  Geophys  Res, 
70:3687-3700 


19D.  Krankowsky,  F.  Arnold,  H.  Weider,  J.  Kissel,  and  J.  Zahringer, 
1972,  "Positive  Ion  Composition  in  the  Lower  Ionosphere,"  Radio  Sci , 
7:93-98 


20R.  S.  Narcisi,  A.  D.  Bailey,  L.  E.  Wlodyka,  and  C.  R.  Philbrick,  1972, 
"Ion  Composition  Measurements  in  the  Lower  Ionosphere  During  the  November 
1966  and  March  1970  Solar  Eclipses,"  J  Atmos  Terr  Phys,  34:647-658 


21A.  Johannessen  and  D.  Krankowsky,  1974,  "Daytime  Positive  Ion  Composition 
Measurement  in  the  Altitude  Range  73  -  137  km  above  Sardinia,"  J  Atmos 
Terr  Phys,  36:1233-1247 

22F.  Arnold  and  D.  Krankowsky,  "A  New  Concept  for  the  D  Region  Ion  Chemistry 
as  Inferred  from  a  Mass  Spectrometer  Measurement,"  paper  presented  at 
International  Symposium  on  Solar-Terrestrial  Physics,  COSPAR,  Sao  Paulo, 
Brazil,  1974 
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Figure  3.  Computed  positive  ion  density  profiles  for  species  that  have  NO 
as  their  core  ion.  The  total  positive  ion  profile  is  given  by 
the  heavy  curve  marked  "T." 
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Figure  4.  Computed  positive  ion  density  profiles  for  species  that  have 

as  their  core  Ion.  The  total  positive  Ion  profile  Is  given  by 
the  heavy  curve  marked  "T." 
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One  point  to  be  noted  is  that  the  code  predicts  H^0+  (H^O)  to  be  the 

dominant  hydronium  ion  in  the  80  to  90  km  region.  This  prediction  is 
because  of  the  proposed  reaction 

N0+  (H20)2  +  hv  +  H30+  +  HN02 

and  the  subsequent  hydration  of  the  H^O  .  Attempted  measurements23  of 

the  photon  absorption  cross  section  of  N0+  (H„0)9  have  indicated  that  it 

-20  2  c  ^ 

is  below  the  detection  threshold  of  10  cm  .  However,  our  modeling  cal¬ 
culations  show  that  the  above  reaction  forms  a  link  between  N0+  and  H90+ 

-212  J 

even  when  a  cross  section  as  small  as  2  x  10  cm  is  assumed.  Specifi¬ 
cally,  at  80  km  the  ratio  of  H30+  (H20)2/H30+  (H20)  as  computed  from  the 

model  is  0.06,  1.3,  and  12  for  photodestruction  cross  sections  of  2  x 
-19  2  -21  2  2 

10  cm  ,  2  x  10  cm  ,  and  0.0  cm  ,  respectively.  Therefore,  relative 

ratio  of  the  hydrates  of  H30+  is  a  sensitive  indicator  of  the  importance 

of  the  photolysis  of  N0+  (H20)2  (whereas  the  total  positive  ion  density 

varies  only  a  few  percent  for  the  changes  in  cross  section  used  above). 
Rocket  ion  measurements  are  ambiguous  at  this  point.  Some  have  shown  an 

abundance  of  H30+  and  H30+  (H20)  ions  in  this  altitude  range,11*24  while 

others  have  indicated  that  H30+  (H20)2  is  the  dominant  hydronium  ion 

cluster.25  The  major  experimental  difficulty  lies  in  the  determination  of 
the  degree  of  breakup  for  the  cluster  ions.  Resolution  of  this  experi¬ 
mental  problem  could  yield  an  answer  to  the  question  of  whether  the 

photolysis  of  N0+  ( H20 ) 2  actually  occurs  in  the  D  region. 


23G.  P.  Smith,  P.  C.  Cosby,  and  J.  T.  Moseley,  1977,  "Photodissociation 
of  Atmospheric  Positive  Ions,  I;  5300-6700A,  •  ■  Chem  Phys,  67:3818-3828 

nR.  S.  Narcisi  and  A.  D.  Bailey,  1965,  "Mass  Spectrometric  Measurements 
of  Positive  Ions  at  Altitudes  from  64  to  112  Kilometers,"  J  Geophys  Res, 
70:3687-3700 

24P.  A.  Zbinden,  M.  A.  Hidalgo,  P.  Eberhardt,  and  J.  Geiss,  1975,  "Mass 
Spectrometer  Measurements  of  the  Positive  Ion  Composition  in  the  D  and 
E  Regions  of  the  Ionosphere,"  Planet  Space  Sci,  23:1621-1642 

25F.  Arnold  and  D.  Krankowsky,  1977,  "Ion  Composition  and  Electron  and 
Ion  Loss  Processes  in  the  Earth's  Atmosphere,"  Dynamical  and  Chemical 
Coup! ing,  93-127,  B.  Grandal  and  J.  Holtet,  eds.,  Dr.  ReideT,  Dordrecht 
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Figure  5  illustrates  the  summed  profiles  for  the  two  main  ion  groupings-- 

the  one  with  N0+  as  the  core  ion  and  the  other  with  the  as  the  core 

ion.  The  crossover  point  between  the  two  groupings  is  75  km,  which  is  in 
reasonable  agreement  with  experimental  observations.  Figure  5  shows  that 
the  region  of  60  to  75  km,  where  the  simulated  total  positive  ion  profile 
is  less  than  the  experimentally  determined  one  (see  figure  2),  is  clearly 

dominated  by  H30+  cluster  ions.  As  mentioned  in  the  previous  section, 

one  possible  reason  for  the  lower  simulated  ion  densities  is  that  the 
effective  recombination  rate  may  be  too  large.  This  large  recombination 
rate  may  in  turn  have  two  causes.  One  cause  may  be  that  the  code  pre¬ 
dicts  multiple-cluster  ions  which  are  larger  (in  mass)  than  the  ones  which 

were  actually  present.  If  the  dominant  cluster  ion  were  H30+  (H20) 

rather  than  the  predicted  H30+  (H20)3,  the  electron  ion  recombination 

coefficient  would  be  half  as  large.  Nominally  this  would  increase  both 
the  ion  and  electron  densities  approximately  40  percent,  which  would 
improve  the  comparison  between  the  ion  profiles  but  increase  the  discrep¬ 
ancy  between  the  electron  profiles.  Other  experimental  determinations  of 
exact  ion  composition  are  still  subject  to  interpretation.  Some  measure¬ 
ments20*24  indicate  that  H30+  (H^O)  may  be  the  predominant  positive  ion 

in  the  lower  D  region,  while  other  experiments11*25  show  that  the  heavier 
clusters  dominate.  Until  the  question  of  cluster  breakup  is  resolved,  we 
adopt  the  position  that  the  code  correctly  predicts  the  heavier  positive 
ion  clusters.  The  other  cause  of  a  too-large  effective  recombination 
rate  may  be  in  the  negative  ion  concentrations  and  negative  ion  to  elec¬ 
tron  ratio.  As  a  "rule-of-thumb,"  the  ion-ion  recombination  coefficients 
are  one  to  two  orders  of  magnitude  smaller  than  electron-ion  recombination 


20R.  S.  Narcisi,  A.  D.  Bailey,  L.  E.  Wlodyka,  and  C.  R.  Philbrick,  1972, 
"Ion  Composition  Measurements  in  the  Lower  Ionosphere  During  the  November 
1966  and  March  1970  Solar  Eclipses,"  J  Atmos  Terr  Phys,  34:647-658 

24P.  A.  Zbinden,  M.  A.  Hidalgo,  P.  Eberhardt,  and  J.  Geiss,  1975,  "Mass 
Spectrometer  Measurements  of  the  Positive  Ion  Composition  in  the  D  and 
E  Regions  of  the  Ionosphere,"  Planet  Space  Sci ,  23:1621-1642 

nR.  S.  Narcisi  and  A.  D.  Bailey,  1965,  "Mass  Spectrometric  Measurement 
of  Positive  Ions  at  Altitudes  from  64  to  112  Kilometers,"  J  Geophys  Res, 
70:3687-3700  - 

25F.  Arnold  and  D.  Krankowsky,  1977,  "Ion  Composition  and  Electron  and 
Ion  Loss  Processes  in  the  Earth's  Atmosphere,"  Dynamical  and  Chemical 
Coup! ing,  93-127,  B.  Grandal  and  J.  Holtet,  eds..  Dr.  Reidel,  Dordrecht 
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Figure  5.  N0+  cluster  ions  and  H^O  cluster  ions  represent  the  sums  of  the 

profiles  of  figures  3  and  4,  respectively.  Below  65  km  the  total 
positive  ions  are  very  well  approximated  by  the  H30+  core  cluster 
ions,  above  85  km  by  tbe  N0+  core  cluster  ions. 
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coefficients.26*27  Thus,  an  appreciable  concentration  of  negative  ions  in 
the  60  and  75  km  region  would  increase  the  simulated  positive  ion  densities 
and  thus  improve  the  agreement  of  the  two  positive  ion  profiles  in  figure 
2,  while  also  possibly  reconciling  some  of  the  differences  between  the 
simulated  and  experimentally  determined  electron  profiles.  This  idea  will 
be  explored  further  in  the  next  section. 


NEGATIVE  ION  CHEMISTRY 

The  major  formation  process  for  the  initial  negative  ion  is  the  three 
body  attachment  of  electrons  to  02 

e  +  02  +  M  -*  0^  +  M. 

However,  there  are  also  rapid  detachment  processes  for  0^ 

02  +  0  -*•  0g  +  e, 

°2  +  02(\)*  2°2  * 


which  tend  to  limit  the  02  density  in  the  upper  D  region.  The  next  impor¬ 
tant  step  in  the  formation  of  more  complex  and  stable  negative  ions  is  a 
charge  exchange  reaction  between  03  and  02 

©2  +  O3  -*■  03  +  02 , 


26D.  Smith,  N.  G.  Adams,  and  M.  J.  Church,  1976,  "Mutual  Neutralization 
Rates  of  Ionospherical ly  Important  Ions,"  Planet  Space  Sci,  24:697-703 

270.  Smith  and  M.  J.  Church,  1977,  "Ion-Ion  Recombination  Rates  in  the 
Earth's  Atmosphere,"  Planet  Space  Sci,  25:433-439 
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followed  by  the  reaction 


0^  +  C02  -*  CO3  +  o2. 

The  COj  ion,  however,  is  also  dissociated  by  atomic  oxygen 
CO"  +  0  -  o2  +  co2 

which  short  circuits  the  negative  ion  formation  process  back  to  the  02 
ion.  The  next  major  step  is  the  reaction 

CO"  +  NO  -*  NO"  +  C02 

which  is  followed  by  the  reaction 

NO"  +  03  -*•  NO"  +  02 

to  form  the  negative  ion  with  the  largest  electron  affinity,  NO^. 

The  major  paths  along  which  the  charge  flows  are  illustrated  in  figure  6. 
Although  there  are  numerous  other  reactions  which  enter  in,  the  main 
sequence  has  been  outlined  above.  Of  interest  is  the  secondary  path  of 

0“  clustering  to  form  0^  which  in  turn  switches  with  C02  to  form  CO^. 

While  this  is  not  a  major  path  for  the  flow  of  charge,  appreciable  buildup 

of  the  CO"  ion  can  occur  because  of  its  small  photodissociation  cross 

section.28*29  Hydration  reactions  occur  for  essentially  all  of  the 


28G.  P.  Smith,  L.  C.  Lee,  P.  C.  Cosby,  J.  R.  Peterson,  and  J.  T.  Moseley, 
1978,  "Photodissociation  and  Photodetachment  of  Molecular  Negative  Ions, 

V:  Atmospheric  Ions  from  7000  to  8400  A,"  J  Chem  Phys,  68:3818-3822 

29M.  L.  Vestal  and  G.  H.  Mauclaire,  1977,  "Photodissociation  of  Negative 
Ions  Formed  in  C02  and  C02/02  Mixtures,"  J  Chem  Phys,  67:3758-3766 
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NEGATIVE  ION  FLOW  CHART 

Figure  6.  Major  paths  for  flow  of  negative  charge  down  the  negative  ion 

chain.  Solid  boxes  indicate  ions  with  appreciable  concentrations; 
dashed-line  boxes  indicate  important  intermediary  ions  which  have 
very  small  concentrations. 


negative  species,  with  hydrated  forms  of  the  ions  expected  to  be  dominant 
at  lower  altitudes.  Chemical  considerations  indicated  that  nitric  acid, 

HNOg,  will  replace  water  in  the  NOj  core  ion  clusters  in  the  strato¬ 
sphere.30’31 

The  negative  ion  chemistry  as  it  is  currently  understood  may  be  divided 
into  conveniently  separable  phases.  The  initial  entry  into  the  negative 
ion  chain  is  through  the  three-body  attachment  of  electrons  to  02>  The 

major  negative  ions,  which  also  serve  as  core  ions  for  clusters,  are  02> 

COj,  and  NOj.  Atomic  oxygen  is  the  major  inhibitor  of  charge  flow  down 

the  chain  because  it  attacks  both  02  and  COj.  Ozone  is  the  major  expeditor 

of  flow  down  the  chain  by  bridging  the  gaps  from  02  to  CO^  and  then,  with 

the  aid  of  NO,  to  NO^.  Clustering  becomes  important  as  one  descends  in 

altitude  and  cluster  ions  will  be  the  dominant  form  of  ion  in  the  lower 
mesosphere  and  throughout  the  stratosphere. 

The  computer  code  simulations  of  the  negative  ion  densities  are  illustrated 
in  the  set  of  figures  7  through  10.  Figure  7  shows  the  altitude  profiles 

for  those  ions  which  have  02  as  their  core  ion.  The  figure  also  shows  the 
0"  ion  density  profile  since  a  major  source  of  0"  is  charge  exchange  with 
Og.  0"  and  02  are  predicted  to  be  the  major  negative  ions  in  the  upper  D 
region,  although  in  terms  of  the  total  charge  they  are  minor  negative  spe¬ 
cies.  C0^  is  considered  to  be  an  02  core  ion  cluster;  and  CO^  and  its 
water  cluster  follow  the  02  profile  very  well  in  the  65  to  40  km  region. 
Below  40  km,  three-body  clustering  reactions  become  so  efficient  that  CO^ 
and  its  hydrate  actually  begin  to  increase  at  the  expense  of  the  simpler 
02  ion.  Figure  8  shows  the  altitude  profiles  of  the  CO3  ion  and  its  water 
cluster.  The  CO^  density  once  again  follows  the  02  rather  closely,  but 
decreases  less  rapidly  than  the  02  density  in  the  stratosphere  due  to 
the  increasing  ratio  of  ozone  as  one  descends  in  altitude  in  this  region. 
CO'  is  computed  to  be  the  most  abundant  negative  ion  in  the  60  and  65  km 
region.  NOj  rapidly  becomes  the  dominant  negative  ion  below  60  km  as 
is  shown  in  figure  9,  with  its  clusters  of  water  and  nitric  acid  becoming 


30F.  C.  Fehsenfeld,  C.  J.  Howard,  and  A.  J.  Schmeltekopf ,  1973,  "Gas  Phase 
Ion  Chemistry  of  HNO3,"  J  Chem  Phys,  63:2835-2841 

31J.  C.  Harris,  D.  G.  Moos,  N.  R.  W.  Swann,  G.  F.  Neill,  and  P.  Glldwarg, 
1976,  "Simultaneous  Measurements  of  H20,  N02,  and  HN03  in  the  Daytime 

Stratosphere  from  15  to  35  km,"  Nature,  259:300-301 
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Figure  7.  Profiles  of  negative  ions  that  have  0^  as  their  core;  C0^  is 

taken  to  be  O2  (CO2).  0  profile  is  also  shown.  The  total 

positive  ion  profile  is  given  by  the  heavy  line  marked  "T." 

Above  55  km  the  total  negative  ion  profile  is  shown  by 

Below  55  km  the  total  positive  and  total  negative  profiles  overlap 
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Figure  8.  Profiles  of  the  negative  ions  that  have  COj  as  their  core  ion. 

The  total  negative  Ion  profile  is  shown  by  the  heavy  line  marked 
above  55  km  and  "T"  below  55  km. 
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Figure  10. 
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.  Profiles  of  the  negative  ions  that  have  NO3  as  their  core  ion. 

^gative  Ion  profile JS  shown  by  the  heavy  line  marked 
above  55  km  and  "T "  below  55  km. 


0,  C,  and  N  represent  the  profiles  of  the  sums  of  the  negative 
ions  that  are  given  in  figures  7,  8,  and  9,  respectively.  The 
electron  density  profile,  E,  is  also  shown. 
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the  more  important  ions  in  the  stratosphere.  Only  one  molecule  each  of 

water  and  nitric  acid  were  clustered  to  NO^  in  this  model.  More  complex 

and  mixed  clusters  have  not  been  included  at  the  present  but  most  probably 
do  exist  in  nature. 

Figure  10  shows  the  three  groups  of  negative  ions,  0^,  CO3  and  NO^,  and 
the  total  charge.  These  groupings  do  not  show  the  sharp  distinctions 
with  altitude  that  the  N0+  and  HjO+  groups  did  in  figure  5,  but  the  tran¬ 
sitions  from  0^  to  CO"  to  NO^  groupings  are  still  apparent  as  one  descends 

in  altitude.  Table  1  lists  the  dominant  positive  and  negative  species  and 
their  corresponding  altitude  ranges.  Certain  features  should  be  noted  when 
figures  5  and  10  are  compared.  First,  clustering  to  positive  ions  is  an 
important  process  throughout  the  entire  D  region,  whereas  clustering  to 
negative  ions  only  becomes  dominant  below  60  km.  Secondly,  the  total 
negative  ion  density  quickly  decreases  above  60  km.  These  features  are 

due  mainly  to  the  fact  that  the  negative  ions  0^  and  CO3  are  attacked  by 

atomic  oxygen,  whose  density  increases  with  altitude,  and  to  the  fact  that 
the  formation  of  the  more  stable  negative  ions  depends  on  ozone,  whose 
density  decreases  with  altitude.  The  DAIRCHEM  computer  code  therefore 
predicts  small  negative  ion  densities  in  the  D  region. 

Comparison  of  the  code’s  predictions  with  available  experimental  data 
is  somewhat  disquieting.  Mass  spectrometer  measurements  have  indicated 
the  presence  of  negative  ions  throughout  the  D  region,32*33*20  although 
the  absolute  concentrations  are  still  in  question.  Rapid  decreases  in 
the  negative  ion  density  with  increasing  altitude  have  been  observed  in 
the  vicinity  of  78  and  90  km,  whereas  the  code  predicts  that  the  negative 
ion  density  begins  to  decrease  above  the  much  lower  altitude  of  60  km. 

The  presence  of  heavy  negative  ion  clusters,  often  with  masses  greater  than 
150  u,  has  been  experimentally  inferred  in  the  70  to  90  km  region.20*32 

The  code  predicts  only  the  simple  negative  ions,  0^,  0",  CO3,  and  N0j, 


32R.  S.  Narcisi,  A.  D.  Bailey,  L.  Della  Lucca,  C.  Sherman,  and  D.  M.  Thomas, 
1971,  "Mass  Spectrometric  Measurements  of  Negative  Ions  in  the  D  and  Lower 
E  Regions,"  J  Atmos  Terr  Phys,  33:1147-1159 

33F.  Arnold,  J.  Kissle,  D.  Krankowsky,  H.  Wieder,  and  J.  Zahringer,  1971, 
"Negative  Ions  in  the  Lower  Ionosphere:  A  Mass  Spectrometric  Measurement," 

J  Atmos  Terr  Phys,  33:1169-1175 

20R.  S.  Narcisi,  A.  D.  Bailey,  L.  E.  Wlodyka,  and  C.  R.  Philbrick,  1972, 

"Ion  Composition  Measurements  in  the  Lower  Ionosphere  During  the  November 
1966  and  March  1970  Solar  Eclipses,"  J  Atmos  Terr  Phys,  34:647-658 
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TABLE  1.  QUIET  TIME  D  REGION 


in  this  region.  The  electron  and  total  positive  ion  density  profiles  shown 
in  figure  2,  which  were  derived  from  independent  rocket  measurements,  imply 
an  increasing  negative  ion  density  in  the  60  and  85  km  region  which  equals 
and  even  exceeds  the  electron  density.  Again  the  known  ion  chemistry  pre¬ 
dicts  a  decreasing  negative  ion  density  in  this  region. 

It  has  become  evident  at  this  point  that,  in  contrast  to  the  positive  ion 
chemistry,  there  are  several  major  uncertainties  in  connection  with  the 
negative  ion  chemistry  and  the  resultant  electron  and  negative  ion  den¬ 
sities.  Several  suggestions  can  be  made  which  may  help  to  reconcile  or 
at  least  better  define  the  current  uncertainties. 

When  the  rocket  measurements  yielding  the  data  shown  in  figure  2  were  made, 
there  were  no  simultaneous  measurements  of  ozone  and  atomic  oxygen  made. 
Considering  the  importance  of  these  minor  neutral  species  in  the  formation 
of  negative  ions,  detailed  profiles  in  the  60  to  90  km  region  would  be 
helpful.  The  DAIRCHEM  code  does  determine  representative  [0^]  and  [0] 

profiles  based  on  time-dependent  photochemical  considerations,  so  that  the 
major  effects  are  already  accounted  for.  Detailed  comparisons  with  experi¬ 
mental  profiles  would  help  to  resolve  questions  concerning  the  exact  alti¬ 
tudes  of  negative  ion  formation,  but  would  not  resolve  major  discrepancies 
involving  total  negative  ion  concentrations. 

The  comparison  of  computer  code  results  and  experimental  data  would  strongly 
imply  that  there  is  another  method  of  negative  ion  formation  other  than  the 

known  three-body  attachment  of  electrons  to  Og.  The  difference  between  the 

measured  electron  and  total  positive  ion  densities  in  figure  2  also  suggests 
that  a  negative  species  not  currently  modeled  is  present  in  the  D  region. 

The  presence  of  ice  particulates  of  submicron  size  has  been  suggested34 
as  a  means  of  providing  an  alternate  electron  attachment  mechanism  and 
negative  species.  Analysis35  and  modeling36  of  data  from  the  1966  solar 
eclipse  have  yielded  similar  conclusions. 37  The  previously  mentioned 


34E.  T.  Chesworth  and  L.  C.  Hale,  "Ice  Particulates  in  the  Mesosphere," 
1974,  Geophys  Res  Letters,  1:347-350 

35R.  D.  Sears,  1972,  "Analysis  of  the  1966  Solar  Eclipse  Data,"  LMSC 
D246526  (DNA  2863F),  Lockheed,  Palo  Alto  Research  Laboratory 

36M.  G.  Heaps,  F.  E.  Niles,  and  R.  D.  Sears,  1978,  "Modeling  the  Ion 
Chemistry  of  the  D  Region:  A  Case  Study  Based  Upon  the  1966  Total  Solar 
Eclipse,"  ASL-TR-0015,  Atmospheric  Sciences  Laboratory,  White  Sands  Missile 
Range,  NM 

37M.  G.  Heaps,  "Inclusion  of  Small  Particulates  to  Explain  Variations  of 
D  Region  Electron  Densities  Observed  During  a  Solar  Eclipse,"  ASL-TR 
(In  preparation) 


measurements  of  heavy  negative  ions  also  lend  indirect  support  to  the 
possible  role  of  small  particulates  in  D  region  processes.  An  area  of 
further  study  would  be  electron  attachment/detachment  processes,  ion 
clustering  and  growth  mechanisms,  and  the  influence  of  nongas-phase 
chemistry  in  the  lower  ionosphere. 


CONCLUSION 

The  positive  ion  chemistry  portion  of  the  DAIRCHEM  computer  code  may  be 
considered  to  be  validated  by  this  comparison.  In  need  of  further  study 
are  electron  attachment/detachment  processes  and  the  possible  role  of 
particulates  in  D  region  ion  chemistry. 
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